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ABSTRACT 


A continuum  theory  of  viscoelasticity  is  developed  which  allows 
non-affine  deformation,  defined  in  an  appropriate  manner.  The 
constitutive  equation  is  a generalization  of  that  obtained  from 
molecular  theory  with  the  addition  of  one  scalar  parameter  which 
becomes  important  for  large  deformations.  The  theory  is  applied  to 
simpler  shear  flows,  the  scalar  parameter  being  determined  to  match 
certain  experimental  data.  The  theory  shows  good  agreement  with 
all  data  examined.  The  paper  concludes  with  the  development  of  a 


general  ncn-affine  thermodynamic  theory. 
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1.  Introduction.  A constitutive  equation  for  viscoelastic  fluids  can 
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be  derived  systematically  from  the  molecular  theory  of  Gaussian  networks.  The 
same  constitutive  equation  can  be  derived  from  the  molecular  bead-spring  model 
with  Hookean  springs.  (See  [1]  for  references).  It,  hence,  has  a sound 
theoretical  foundation.  Unfortunately,  it  does  not  describe,  even  qualita- 
tively, some  significant  viscoelastic  behavior  such  as  the  '‘snear  thinning" 
of  viscosity. 

One  of  the  assumptions  of  the  molecular  theory  of  networks  is  that  the 
junctions  move  affinely  with  an  equivalent  continuum*  The  anologous  assump- 
tion of  the  bead-spring  model  is  that  the  surrounding  Newtonian  fluid  moves 
affinely  with  an  equivalent  continuum.  It  is  our  object  to  relax  the 
affine  assumption  present  in  both  models  by  developing  a continuum  theory 
that  allows  two  histories  of  deformation  which  may  be  non-affine.  One 
history  is  made  to  give  rise  to  the  current  state  of  stress  according  to  the 
constitutive  equation  from  molecular  theory.  The  second  is  made  the  observed 
smooth  continuum  deformation.  The  two  motions  are  connected  by  an  appro- 
priate constitutive  equation. 

In  the  next  section  we  develop  the  constitutive  equation  which  allows 
non-affine  deformation.  It  is  just  as  easy  to  use  as  that  based  on  the 
affine  assumption,  but  is  able  to  describe  viscoelastic  behavior  consider- 
ably better  for  all  the  data  we  have  examined.  This  improvement  is 
purchased  by  the  addition  of  one  scalar  parameter.  The  constitutive  equa- 
tion is  then  used  to  examine  shear  flow  and  the  new  scalar  parameter  in 
the  theory  determined  to  match  certain  experimental  data.  Steady,  sudden 
start  up  and  small  oscillatory  shear  flows  are  examined.  It  is  not  our 
intention  here  to  make  an  intensive  comparison  with  data  or  with  the  results 
of  the  many  empirical  models  found  in  the  literature.  Such  comparisons 

•Actually  it  is  sufficient  to  assume  only  that  boundary  points  of  the  network 
move  affinely.  It  can  then  be  proved  that  ensemble  - average  positions  of 
all  junctions  move  affinely  [1]. 
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will  follow  In  a later  paper.  We  here  wish  to  develop  a sound  theoretical 
basis  for  the  theory  and  to  compare  results  with  enough  data  to  indicate  its 
potential  usefulness.  In  this  regard  we  develop  in  the  last  section  of  this 
paper  a thermo-viscoelastic  theory  within  the  framework  of  generalized 


continuum  mechanics.  This  theory  contains  the  theory  of  the  next  section  as 
a special  case. 


2.  Development  of  the  constitutive  equation.  Consider  two  deformation 
histories  which  we  denote  by  a and  0.  The  locations  of  particles  in  these 
histories  are  described  by  Cartesian  coordinates  x^  and  respectively. 

If  the  two  histories  are  related  by  a compatible  deformation,  there  is  a 
one-to-one  correspondence 


<xr  x2,  x3,  t) 


(2.1) 


describing  this  deformation  at  every  time  t.  The  local  deformation  is 
described  by 


d5i  ■ Tu 


(2.2) 


where 


rij  ■ V9’j 


(2.3) 


are  the  deformation  gradients.  Define  displacement  component's  r^  by 


5i  ■ ‘i  + ri 


(2.4) 


The  velocity  components  of  particles  in  each  history,  x^  and  are  related  by 


‘ 4i  + *i 


(2.5) 


Differentiation  of  (2.5)  gives 

I 

A F ■ x + r 
lr  rk  l,k  l,k 


(2.6) 


where 


AU  " 


(2.7) 


and  the  comma  denotes  differentiation  with  respect  to  x^. 

We  next  derive  transformation  rules  under  a change  of  frame  of  reference 
for  the  various  kinematic  quantities  defined  above.  The  relative  orientation 
of  the  two  frames  can  be  described  by  a translation  c(t)  and  an  orthogonal 
rotation  tensor  Q(t).  Let  be  Cartesian  coordinates  associated  with 

the  second  frame  and  and  be  components  taken  with  respect  to  base 

vectors  in  the  first  frame.  It  can  be  shown  that 

xi  * Vj  + Ci 

Differentiation  of  this  relation  gives: 

’ Vj  + Vj  + 6i 

Since  Q is  orthogonal, 

^ijQik  + Qij^ik  ’ ° 


If  the  two  frames  coincide  at  current  time,  = 6 and  the  above  give: 


Xi  " Xi  + Ci 


(2.8) 


xi  " xi  + Qijxj  + Ci 


(2.9) 


QU  + QJi 


(2.10) 


and  also  transform  according  to  (2.8)  and  (2.9)  so  that  the  rules 

for  r^  and  r^  are: 


" 'i  ’ W Vd 


(2.11) 


Differentiation  of  (2.9)  and  (2.11)  give  the  rules: 


Xi,k  " xi,k  + Qik 


*i,J  " ri,j  + ^ikFkj  ' ^ij 


(2.12) 


(2.13) 


From  (2.6),  (2.12)  and  (2.13)  we  obtain 


Aij  " Aij  + $ij 


(2.14) 


If  histories  a and  B are  identical,  we  say  they  are  affine.  Then, 

I 

P . ■ 6 , , r,  “0  and  A,,  * x , . Non-affine  histories  may  be  incompatible 
lj  lj  A 1 * J 

in  which  case  a global  correspondence  (2.1)  is  not  defined.  Eqs.  (2.3), 
(2.4),  (2.5)  and  (2.7)  then  cease  to  have  meaning.  Ue  still  make  use  of 
fields  r^,  x^  and  which  satisfy  Eq.  (2.6)  and  transformation 

relations  (2.8)  through  (2.14).  History  a,  which  is  assumed  to  be  smooth 
and  homogeneous,  is  taken  as  the  history  observed  on  the  macroscopic  scale. 
History  B will  be  taken  to  give  rise  to  the  stress  field  in  a manner  to 
be  made  precise  presently.  F represents  the  local  deformation  of  a to 
B according  to  Eq.  (2.2)  and  A represents  the  velocity  gradient  of  history 
B in  a local  sense  also. 

If  tensor  A is  made  to  depend  on  x.  , by  making  F and  r,  , 
depend  on  x^  j , there  will  result  a continuum  model  which  represents  a 
simple  fluid.  We  develop  a theory  of  this  type  based  on  the  following 
assumptions : 

(a)  Fij(t)  = 6 . At  current  time  the  spatial  labelling  of  both 

histories  is  the  same. 

(b)  r,  is  a linear  isotropic  function  of  x . 

(c)  The  Cauchy  stress  T depends  on  the  history  of  tensor  A. 

(d)  The  deformation  is  incompressible  according  to  the  equation: 

ilfl  “ 0 (2.15) 

By  assumptions  (a)  and  (b) , the  velocity  gradient  A becomes  a function 
of  the  gradient  x^  ^ of  a.  Assumption  (c)  says  that  the  stress  depends 
on  the  history  a,  which  is  the  constitutive  relation  of  a simple  fluid. 
Assumptions  (b)  and  (d)  result  in  the  following  form: 


where  b and  c are  constants.  Use  of  assumption  (a)  and  Eqs.  (2.6)  and 


(2.16)  yield: 


Au  ■ (b+1)ii,j  + c *i,i 


(2.17) 


We  require  (2.17),  which  plays  the  role  of  a constitutive  equation,  to  be 
form  invariant  under  a change  of  frame  of  reference.  Substitution  of  (2.12) 
and  (2.14)  in  (2.17)  and  use  of  (2.10)  results  in: 

(c-b)  Q1;j  - 0 

which  must  hold  for  arbitrary  antisymmetric  Q.  Hence,  c * b and  we  write 
both  in  terms  of  another  constant  a,  as  follows: 

_ . h _ a-1 

c b ~T 

Eqs.  (2.16)  and  (2.17)  become 


*i.j  ' (a‘l)  ‘(i.j) 


a+1  . a-1  . 

Aij  " 2 Xi,j  2 xj,i 


(2.18) 


where  x..  is  the  symmetric  part  of  x . . . 

Eq.  (2.18)  says  that  the  antisymmetric  parts  of  A.,  and  x.  . are  the 

1 * J 

same.  That  is,  a and  8 have  the  same  vorticity.  The  symmetric  parts 


are  related  by: 


A(i,j)  " 3 X(i* j) 


(2.19) 


Hence,  the  rate-of-deformation  of  elements  in  8 is  a times  the  rate-of- 


deformation  in  a.  a = 1 corresponds  to  affine  motion,  the  deviation  from 

this  value  measuring  the  deviation  from  affine  motion. 

There  are  various  ways  assumption  (2.16)  can  be  generalized  and  still 

result  in  a simple  fluid  theory,  r can  be  allowed  to  be  a nonlinear 

J 

isotropic  or  anisotropic  function  of  x.  . or  it  can  be  allowed  to  depend 

*■»  J 

on  the  history  of  x . . 
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Assumption  (c)  is  written  in  the  form: 


T + pi  - F {E(t, t) } (2.20) 

X»— co 

where  p is  a scalar  pressure  and  E is  a deformation  measure  defined  as 
the  solution  of : 


9E(t, t 


A(t)  E(t, t) , 


E(T,T)  - I 


(2.21) 


where  9E/9t  is  a material  time  derivative.  It  can  be  shown  that  E also 
satisfies: 


= -E(t  , t)  A(t)  , E(t,  t)  - I 


(2.22) 


With  the  generality  of  Eq.  (2.20)  this  theory  is  equivalent  to  the 
usual  simple  fluid  theory.  Something  new  is  obtained  when  a special  con- 
stitutive equation  is  used.  A constitutive  equation  derived  from  molecular 
theory  is  given  by  [1J: 


T + pi  - f m(t-T)K(T,t)K(T,t)T  dt 


(2.23) 


where  K is  the  gradient  9x(t)/9x(x)  and  m a scalar  kernel  func- 

tion. Eq.  (2.23)  is  derived  from  a molecular  network  model  of  Gaussian  chains 
where  the  junctions  are  assumed  to  move  affinely  with  an  equivalent  continuum. 
It  can  also  be  derived  from  a bead-spring  molecular  model  where  the  springs 
are  Hookean  and  the  bead-spring  system  is  assumed  to  be  surrounded  by  a 
Newtonian  fluid  which  moves  affinely  with  an  equivalent  continuum.  One  of 
the  most  important  features  of  this  model  is  that  it  predicts  a constant 
viscosity  in  steady  sheaf  flow,  a result  in  disagreement  with  most  measure- 
ments. We  note  that  tensor  K satisfies  the  equations  (2.21)  and  (2.22) 

with  E replaced  by  K and  A replaced  by  velocity  gradient  i.  . of 

1*  J 

deformation  a. 


1 

! 

i 


! 


A continuum  model  which  relaxes  the  affine  assumption  is  achieved  by 
identifying  the  equivalent  continuum  of  the  models  with  the  history  a but  allow- 
ing the  molecular  models  to  move  affinely  with  history  3 by  replacing 
K in  Eq.  (2.23)  with  E. 

t + 

T + pi  - / m(t-T)E(T,t)E(T,t)T  dT  (2.24) 

—OO 

In  what  follows  the  consequences  of  (2.24)  will  be  investigated.  There  is 
one  additional  parameter,  a,  in  this  new  theory  as  compared  with  the  affine 
molecular  models  which  lead  to  (2.23).  The  choice  of  this  parameter  allows 
better  agreement  with  experimental  results  than  can  be  obtained  with  the 
affine  model.  In  particular,  the  viscosity  is  now  found  to  vary  realisti- 
cally with  the  shear  rate. 

An  alternate  constitutive  equation  is  obtained  by  integrating  (2.24) 
by  parts,  to  obtain 

T + pi  - / G(t-T)E(T,t)D(T)E(x,t)T  dT  (2.25) 


T 


where 

m(t)  - - \ ~ G(t) 

and 

D - | (A  + A+)  (2.26) 

is  the  rate-of-deformation  tensor  of  6.  Note  that  D * a x.  . from 
Eq.  (2.19).  When  a =■  1 the  new  model  agrees,  of  course,  with  the  affine 
molecular  model  (2.23).  As  a -*■  0,  it  becomes  the  co-rotational  model  of 
Goddard  and  Miller  [7]*that  is  discussed  by  Bird,  Hassager  and  Abdel  - Khalik 
[3].  In  what  follows  we  assume  that  0 < a < 1 because  it  is  found  that 
values  of  a in  this  range  yield  good  results  when  compared  with  experi- 

ment  whereas  negative  values  of  a do  not  yield  good  results.  Note  that 

*In  taking  the  limit,  Ga  is  held  fixed  as  the  kernel  in  the  co-rotational 
model. 
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a • -1  corresponds  to  a so  called  covariant  model  (equation  2 of  table  6.4, 


ref.  [1]). 


3.  Simple  shear  flow.  We  next  examine  the  results  predicted  by  the 


non-affine  model  for  simple  shear  flow.  This  flow  is  defined  by: 
Xj^  = <(t)  x2  , x2  = x3  * 0 

the  velocity  gradient  and  rate-of-deformation  tensors  are 


(3.1) 


(3.2) 


0 10 


D«-a  K 1 0 0 


0 0 0 


(3.3) 


Integrating  system  (2.21)  yields 


E = E 
11  22 


cos(A  s) 


*12  " 8in(* 


sin  (A  s) 


E31  “ E13  = E32  “ E23  “ 0 ’ E33  " 1 


(3.4) 


where 


1 / 2 ' ^ 

A - ± / l-aZ  ,s  - / (c(t')  df 
2 T 


(3.5) 


The  constitutive  Eq.  (2.25)  yields: 


-8- 


For  the  sudden  start  of  shear  defined  by  <(t)  = 0 
<(t)  ■ k - constant  for  t > 0,  we  obtain  for  t > 0, 


for  t < 0 and 


t 

/ G(t)  cos(2Akx)  dr 
0 


T 


11 


t22  - fx  f G(T)  s±n(2AkT)  dT 


1 


(3.8) 


Note  that  the  first  maximum  of  occurs  at  strain  kt  = tt/4A  and  is 

Independent  of  shear  rate.  In  this  theory  the  parameter  a can  be  chosen 
to  cause  the  first  maximum  to  occur  at  any  given  strain,  whereas  in  the 
co-rotational  limit  the  first  maximum  must  occur  at  kt  = tt/2 . The  affine 
molecular  model,  with  A =»  0,  does  not  predict  overshoot.  Examination  of 
same  data  for  solutions  given  in  [3]  indicates  that  the  maximum  occurs  at 
a strain  of  2.9  to  3.4,  which  corresponds  to  A between  0.23  and  0.27  and 
a between  0.84  and  0.89.  For  a melt  data  [8]  indicates  a maximum  at  a 
strain  of  4.7  which  corresponds  to  A = 0.17  and  a = 0.94.  We  note  that 
our  model,  as  does  the  co-rotational  model,  predicts  stress  components  that 
undergo  an  infinite  number  of  damped  oscillations  with  time.  For  normal 
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stress  differences,  the  first  maximum  occurs  at  twice  the  strain  as  the  first 
maximum  of  shear  stess. 

For  steady  shear  flow,  k * k = constant  for  all  time,  we  obtain  s = 
k(t-x)  and  viscometric  functions: 

T OO 

n(k)  - ~ a f G(t)  cos(2Xkx)  dx 
k z 0 

OO 

V (k)  = --,-22  = -£r  / G(X)  sin(2Xkx)  dx 

1 k2  2Ak  Q 

T -T 

22  33  1 

*2(k)  = 2 ~ = " 2 (1_a)  Vk) 

k 

The  present  model  predicts  negative  which  agrees  with  most  experimental 
results.  If  we  take  a = 0.87  for  solutions  and  a = 0.94  for  melts,  we 
obtain  4^  = -0.07  4^  and  4^  = -0.03  4^,  respectively. 

Any  of  relations  (3.9),  which  are  Fourier  transforms,  can  be  inverted 
to  yield  expressions  for  kernel  G in  terms  of  n,  4^  or  V ^ In 
principle,  we  can  use  these  relations  to  determine  G from  data  on  n»  4^ 
or  ^2’  H the  experimental  material  behaves  exactly  according  to  our  model, 
each  equation  will  yield  the  same  result.  The  co-rotational  model  can  be 
used  in  the  same  manner  to  determine  G.  The  affine  molecular  model  predicts 
constant  n and  4^  and  gives  4^  = 0 for  every  G. 

These  equations  can  be  used  to  derive  relations  between  the  viscosity 
function  n and  the  normal  stress  functions  4^  and  4'2.  For  example, 
from  (3.9) , 


/ G(x) [cos(2Xxx)  - cos(2Xxk)J  dx 
0 


I dx 


0 k2  - x2 


2 a ' 


dx 

f.  .2  2 

0 k - x 
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On  Interchanging  the  order  of  integration  and  using 


( 


f°  cos  (ax)  - cos.(ak)  dx  _ ^r_  sln(ak) 
0 X - TL 

we  obtain 

f°  ^ dx  “ ~ ^(k) 

0 k - x 


(3.10) 


(3.11) 


Eq.  (3.10)  can  be  derived  by  writing  the  integral  as  a Cauchy  principal 
value  and  evaluating  it  using  contour  integration.  Eq.  (3.11)  can  be  made 
to  match  experimental  results  for  polymer  melts  and  for  polymer  solutions 
given  in  [3]  by  taking  X ■ 1/6  and  X - 1/4  respectively*  These  values 
of  X agree  well  with  those  obtained  from  the  shear  start  up  experiments. 

Note  that  the  value  of  the  first  maximum  of  shear  stress  in  sudden  start 
up  is  related  to  the  steady  shear  viscosity  given  by  (3.9)  by 
_ °°  cos(^-  x) 

T - - k / n(kx)  dx  (3.12) 

max  tt  „ , _ 2 


where  k is  taken  to  be  the  same  in  each  flow.  This  equation  is  obtained 
by  inverting  (3.9)  to  obtain  G(t)  in  terms  of  n(k)  and  substituting 
the  result  in  (3.8)^  with  t ■ Tr/4Xk.  This  result  is  independent  of 
material  quantities  G(t)  and  a.  In  particular  this  prediction  agrees 
with  that  of  the  Goddard-Miller  model. 

For  oscillatory  shear  flow,  we  set 
<(t)  * k cos(u)t) 

It 

a * — (8 incot  sincox) 

in  Eqs.  (3.6)  and  (3.7).  If  the  factors  cos(2Xs)  and  sin(2Xs)  are 

expanded  in  series  of  powers  of  k and  the  leading  terms  retained,  we 

obtain  the  following  results  for  small  amplitude  oscillatory  shear  flow: 

*Eq.  (3.11)  was  derived  in  [3]  for  the  co-rotational  model,  X - 1/2,  and 
empirically  modified  to  the  form  (3.11)  in  order  to  fit  experimental  data. 
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' , I 

A'l  * 


r 

r 


T^2  ■ k (n'  cosut  + n"  sinwt) 

T,  - - T . = k^(0  + 6'  cos2ut  + 0"  sin2u>t) 
11  Li  m 


(3.13) 


where 


n'(w)  ■ — a / G(t)  cosut  dx 
1 0 


n"(u>)  » r a / G(x)  sinux  dx 
0 


0 - / G(x)  sinwx  dx 

“ 2U  0 


0'  - y-  f G(x)  (sin2(i)X  - sinwx)  dx 
“ 0 


0"  ■ I G(x)(coswx  - cos2(JX)  dx 
0 


(3.14) 


Comparison  with  Eqs.  (3.9)  shows  that  the  large  deformation  viscosity  functional 
forms  are  related  to  the  oscillatory  forms  as  follows: 


n'(x)  = n(^0 

(3.15) 

n"(x)  = \ x 4'1(^) 

(3.16) 

en(x)  - \ n" (x)  =•  \ Vf) 

(3.17) 

9’(x)  = i In" (2x)  - n"(x)] 

• Vf>  - 2 Vt> 

(3.18) 

e"(x)  » i [n’(x)  - n’(2x)] 

* X - n<x>l 

(3.19) 

Prom  data  on  n'  and  n parameter  X can  be  determined  by  shifting  curve 
n'(x)  horizontally  until  it  coincides  with  curve  n(x)  and  using  Eq.  (3.15). 
It  can  also  be  determined  by  shifting  curve  — n"  to  coincide  with  and 

using  Eq.  (3.16).  The  horizontal  shifts  are  performed  with  the  data  plotted  on 
semi-log  paper,  of  course.  Data  given  in  [3]  for  a polymer  solution  determine 
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X ■ 1/4  from  both  equations.  From  data  for  a melt  given  in  [3]  both  also 

give  X “ 1/4.  This  is  the  one  melt  in  six  discussed  in  [3]  [or  which 

Eq.  (3.11)  also  gives  X * 1/4.  For  the  other  five  (3.11)  yields  X = 1/6. 

For  a soap  solution,  data  given  in  [4]  gives  X = 1/4  for  both  shifts. 

In  conclusion,  all  data  we  have  examined  give  values  of  X between  1/6 

and  1/4.  We  have  found  that  A = 1/4  for  all  polymer  solutions  and  X * 1/6  for 

2 

all  polymer  melts  save  one.  Corresponding  values  of  a are  3/4  and  8/9. 

Note  that,  from  the  form  of  expressions  (3.14),  the  small  amplitude  oscil- 
latory motion  can  be  used  to  determine  the  product  aG,  but  cannot  determine  a 
separately.  In  this  sense,  the  parameter  a is  relevant  to  finite  amplitude 
motions.  Note  that  ^ aG  is  the  "relaxation  modulus"  of  linear  viscoelasticity  [9]. 

4.  A non-affine  general  thermo-viscoelastic  theory.  We  now  formulate 
a theory  within  the  frame  work  of  generalized  mechanics  [2]  that  permits 
a deformation  which  is  not  affine  with  an  equivalent  continuum. 

The  assumptions  of  §2  are  relaxed  so  that  a theory  of  some  generality  is 
achieved  which  includes  the  previous  results  as  a special  case.  This  theory 
is  useful  in  that  it  shows  the  assumptions  made  in  §2  are  consistent 
within  a more  general  framework;  and  it  provides  a proper  setting  for 
extensions  of  the  special  theory. 

The  governing  equations  are  derived  from  the  first  and  second  laws  of 
thermodynamics  by  application  of  appropriate  invariance  requirements.  This 
method  is  due  to  Rivlin  and  Green.  We  do  not  include  kinetic  energy  or  the 
work  due  to  body  forces  in  this  formulation.  As  a result  both  inertia  and 
body  force  terms  are  left  out  of  the  equations  of  motion.  These  effects 
can  be  added  if  one  wishes. 


The  law  of  energy  balance  and  the  Claus ius-Duhem  inequality  are: 


i 


(4.1) 


/ u p dV  * w -•  f qdA 
R A 

/ A p dV  > - / fdA  (4.2) 

R A ° 

where 

* - f [PjXjl  + dA  (4.3) 

A 

is  the  rate-of-work.  p Is  the  stress  vector  conjugate  to  velocity  x of 
history  a.  r Is  a generalized  velocity  field  and  m the  stress  vector 
conjugate  to  it.  x and  r are  identified  with  the  same  quantities  of  §2 
by  assuming  they  transform  according  to  Eqs.  (2.9),  (2.11),  (2.12)  and 
(2.13).  R is  a neighborhood  in  the  - coordinate  space  with  surface  A.  u 
is  the  internal  energy  density,  p the  mass  density,  q heat  flux,  0 temperature  and 
s the  entropy. 

Tetrahedron  analysis  of  (4.1)  yields: 

pi*i  + mi*i  ' q = (Pij*i  + “ij^i  “ qj)mj 
where  P and  M are  stress  tensors,  3 the  heat  flux  vector  and  n the 
unit  normal.  Tetrahedron  analysis  of  (4.2)  gives 

9 - q^  (4.5) 

Under  a translation  of  frame  given  by  (2.9)  and  (2.11)  with  Q = 0 we 
require  that  (4.4)  remain  Invariant.  This  implies  that 


pi  ’ Vj 


Equations  (4.4),  (4.5)  and  (4.6)  imply 


m * M n 
i ij  j 


(4.6) 


(4.7) 


Use  of  (4.4)  and  the  divergence  theorem  in  Eq.  (4.1)  allows  one  to  write 
the  latter  as  a single  volume  integral  from  which  the  local  form  of  the 
first  law  follows. 
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P“  “ Pij,j*i  + PijXi.J  + Mij,jrl  + " qi,i 


For  (A. 8)  to  remain  invariant  under  a translation,  we  must  have 


P - 0 

ij.j 


(4.8) 


(4.9) 


which  is  an  equation  of  motion  expressing  the  balance  of  linear  momentum. 

With  (4.9)  one  term  can  be  deleted  from  (4.8).  For  (4.8)  to  be  invariant 
under  transformations  (2.9),  (2.11),  (2.12)  and  (2.13)  for  all  Q satisfying 
(2.10),  we  must  have 

p + m r+MF  - M 
ij  ik,krj  “ik*jk  "ij 


“ Pji  + Mjk,kri  + MjkFik  " Mji 
which  is  an  expression  of  the  balance  of  angular  momentum. 
The  local  form  of  law  (4.2)  is 

P08  > - qi,i  +iqi  O.i 


(4.10) 


(4.11) 


Introducing  the  free  energy  h * u - 0s  and  eliminating  q.  . between 

If  1 


(4.8)  and  (4.11)  gives: 


pf.<  - p.e  + 

+ Vl.j  - 1 "i  9-i + » 


(4.12) 


where  the  last  term  has  been  added  with  a pressure  multiplier  p to  account 


for  the  assumed  incompressibility  constraint 

^.i  ■ 0 


(4.13) 


As  thermodynamic  extensive  variables  we  use  0,  and  where 


the  latter  are  defined  as  solutions  of  the  systems 

3t  (t,T)  " xi,kHkj>  * H4J> (t.t)  " 6ij 


(4.14) 
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dH 


(2) 


IS1  <C’T>  “ - 6iJ  (4'15) 

The  free  energy  is  assumed  to  depend  on  the  history  of  these  variables  as 


well  as  explicitly  on  0(t). 


h(t)  - h^  (t,  9(t),  0(T),  H(1)(t,T),  H-(t.T)} 

A chain  rule*  for  nonlinear  mappings  yields: 

t 


,<2), 


h(t)  9t  90(t) 


(4.16) 


(1) 


0(t)  + 6h1(t,0(t),0(T),  H(1) , H(2) ; ^— ) 


6E,(t,0(t)f0(T),  H(1),  H(2) ; (4.17) 


3t 


3h  * 

where  Indicates  partial  differentiation  with  respect  to  explicit  depen- 

dence on  t holding  the  other  arguments  fixed  and  6h^  and  6h^  are 
mappings  which  are  linear  in  their  last  arguments.  The  use  of  this  linearity 
and  Eqs.  (4.14)  and  (4.15)  allows  us  to  conclude  that  6h^  is  linear  in 
^(t)  and  (5^2  linear  in  r^  ^(t)  to  obtain  the  forms: 

,U> 


(t,0(t) , 0(x),  H(1),  H(2);  ) 


- 


6h^(t,0(t),0(T),  H(1),  H(2))  i±  (t) 


(4.18) 


«S2  (t,0(t) , 0(T),  H(1) , H(2) ; ^|^- 


(4.19) 


- «h(2)(t,6f(t),  0(X),  H(1),  H(2))  i.  .(t) 

ij 


Chain  rules  of  this  type  are  discussed  in  [5]. 
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Using  (4.17),  (4.18)  and  (4.19),  the  second  law,  Eq.  (4.12),  takes  the  form: 


p ft  + £ qi  e*i  ^ - p<8  + ai|E>>  + Mijtj 
+ <pu  + p - p «h5i})  * 


ij 


ij  ' i,j 


(4.20) 


+ (M  - p 6h<?>)  r. 


ij  " ij  ' "i.j 

If  the  rates  8,  and  r^  are  Independent,  (4.20)  can  hold  for 

all  values  of  the  extensive  variables  only  If 


. . 3h 
8 1(0 


(4.21) 


"u.j  - 0 


Py  + P «y  * P «h<» 


My  - P «h<f 


(4.22) 

(4.23) 

(4.24) 


The  Cauchy  stress  tensor  is  given  by 


. 4 


r'l 


k ■* 

, 


» 


Tij  ' PiJ  + M«  ■ * pSij  + 0 Shy’  + P ^y* 
The  second  law  (4.20)  reduces  to 


3h  .1  . 

P3t+0qi0»i  - 0 


With  (4.22),  Eq.  (4.10)  takes  the  form 


(4.25) 


(4.26) 


pu  + Vjt  " "u  ’ pn  + V»  - 
« 

wWl*  the  first  law  (4.8)  becomes 


(4.27) 


, ~ & 

. k 

v * 

• k 

ij 

f fc 

I * m. 


* " PiJii,J  + Mijri,J  " qi,i  (4.28) 

When  the  rates  in  (4.20)  are  dependent,  we  postulate  that  Eqs.  (4.21)  through 

(4.28)  still  hold.  These  equations  then  imply  that  the  second  law  in  the  form 
(4.20)  holds  but  are  more  restrictive. 
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The  basic  equations  are  (4.9),  (4.13),  (4.14),  (4.15),  (4.21),  (4.22),  (4.23), 
(4.24),  and  (4.28).  To  these  must  be  added  a vector  equation  of  heat  con- 
duction. The  variables  In  these  equations  are  P^,  M^,  x^,  r^,  q^,  s,  0, 

p,  and  H^.  Note  that  p is  considered  known  and  u In  (4.28)  is 

given  in  terms  of  other  variables  by  u ■ h + 6s.  The  constitutive  equations 
(4.23)  and  (4.24)  must  be  compatible  with  (4.26)  and  (4.27). 

In  the  special  theory  of  §2,  we  have 


Fij  " 6ij  ’ Aij  " *i,j  + *i,j 


If  we  take  the  free  energy  functional  (4.16)  in  the  form 


(4.29) 


h(t)  - ■—  f m(t-T)  tr  (E  E+)  dT 


(4.30) 


and  carry  through  the  steps  leading  to  (4.17),  (4.18)  and  (4.19),  we  get 


H ■ 25?  n<0)'  + zfe  If  <M>  « ff  E+>  * 

#hS>“h«>  “lk'jk  dT 


(4.31) 


(4.32) 


Substitution  of  (4,32)  in  (4.25)  yields  the  constitutive  equation  (2.24). 
Substituting  (4.31)  in  inequality  (4.26)  gives  a restriction  on  the  kernel 
function  m. 

5.  Concluding  Remarks.  Under  assumptions  of  §2  the  history  $ 
coincides  with  the  smooth  deformation  a at  current  time,  but  at  other 
times  8 represents  an  Incompatible  deformation  of  a.  In  this  sense  we 
can  say  that  the  deformation  8 is  inhomogeneous.  This  idea  is  similar  to 
that  used  by  Noll  [2]  [4]  in  developing  a theory  suitable  for  describing 
a continuous  distribution  of  dislocations  where  use  is  made  of  an 


i 


- 


inhomogeneous  reference  state.  Noll's  theory  would  seem  to  be  more  suitable 
for  describing  the  behavior  of  materials  with  a preferred  reference  state 
whereas  the  present  theory  is  more  suited  to  fluids. 

We  note  that  in  the  thermodynamic  theory  developed  by  Coleman  (see  [5] 
for  references)  the  free  energy  is  assumed  to  depend  explicitly  on  the 
current  value  of  a strain  measure  and  to  be  differentiable  with  respect  to 
this  measure.  In  the  present  theory  an  assumption  of  this  kind  is  not 
needed,  the  free  energy  functional  (4.16)  depending  only  on  the  history  of 
strain  measures  and  H^.  It  has  been  pointed  out  by  Astarita  and 

Sartl  [6]  that  this  assumption  made  in  the  Coleman  theory  is  a serious  one. 
Implying  that  the  material  can  support  a jump  in  strain,  a deformation  not 
possible  in  a Newtonian  fluid.  Thus  the  Coleman  model  excludes  Newtonian 
fluid  behavior.  Since  this  assumption  is  not  made  in  our  model,  it  may 
be  more  suitable  for  the  description  of  fluid  behavior,  both  viscous  and 
viscoelastic.  The  reason  the  assumption  is  not  needed  here  is  because  of 
the  use  of  strain  measures  ^ and  defined  by  differential 

equations  (4.14)  and  (4.15). 

Finally  we  wish  to  thank  our  colleagues  in  the  Rheology  Research  Center 
for  their  helpful  discussion  of  this  paper.  We  especially  thank 
Professor  R.  B.  Bird  for  his  detailed  suggestions  and  Professor  Arthur  Lodge 
who  provided  the  motivation  for  this  paper  by  pointing  to  the  need  to  relax 
the  affine  assumption  of  the  molecular  network  theory  (Ref.  [1],  pp.  232-235). 
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